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We have studied tunneling-mediated local density of states (LDOS) of the surface layer of a 
bilayer cuprate, where a Zn impurity is located on the second Cu-O layer. When the tunneling 
strength between two Cu-O layers is larger than a critical value, the LDOS on the site just above 
the Zn impurity first exhibits a resonant peak near the Fermi surface. The larger the tunneling 
strength, the stronger the resonant peak. It is also shown that the height of the resonant peak 
oscillates decreasingly with the distance from the site just above the Zn impurity. The location of 
the resonant peak in the surface LDOS depends on doping, energy gap, and the tunneling strength, 
and has an opposite bias voltage to that on its nearest neighboring sites. The results could be 
tested by the STM experiments and be used to further understand the electronic properties of high 
temperature superconductors. 
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High temperature superconductivity in cuprates has 
been the focus of theoretical and experimental investi- 
gations [1-7]. One of the main features of the cuprate 
superconductors is that their electronic properties de- 
pend strongly on the number of Cu-O planes in a unit 
cell. The bilayer LDA band calculations predicted that 
the c-axis electron hopping has an anisotropic form i k = 
— t z (cosk x — cos k y ) 2 / A, where t z is the tunneling strength 
[8]. Recently angle- resolved photoemission spectroscopy 
(ARPES) [2], one direct probe of the electronic struc- 
ture in the momentum space, confirmed the bilayer split- 
ting in the electronic structure of (heavily) overdoped Bi2 
Sr 2 CaCu 2 08+5 (Bi2212) [9,10]. The experimental value 
of t z is of similar magnitude of the superconducting en- 
ergy gap. 

However, scanning tunneling microscopy (STM) [7], 
one direct probe of local density of states (LDOS) at the 
surface of samples in the real space, has not yet reported 
the distinct feature in bilayer cuprates. This is because 
it is very difficult to figure out the bilayer splitting in 
uniform samples by the STM experiments. In Rcf. [11], 
Pan et al. observed a resonant peak near the Fermi sur- 
face on the Zn impurity located at the first Cu-O layer 
of Bi2212 and discovered a four-fold symmetric LDOS 
pattern by the STM experiment. The LDOS pattern 
near the Zn impurity is consistent with the symmetry 
of the superconducting order parameter. These experi- 
mental phenomena were explained by a unitary potential 
and the blocking effect due to the Bi atom just above 
the Zn impurity [12] (For an alternative model, see Ref. 
[13]). In this work, we investigate the surface LDOS in 
bilayer Bi2212, where a Zn impurity is located on the 
second Cu-O layer so that we can find out the features 
of the tunneling-mediated impurity states and the tun- 
neling splitting, which could be observed by the STM 
experiments. 

The Hamiltonian describing the scattering of quasipar- 



ticles from a single impurity and the inter-layer tunneling 
in a d-wave bilayer superconductor can be written as 

H=^2 £ kC+ k(T c„, k(T + ^ A k (c+ kT c+_ ki + c n ,_ ki c n , kT ) 
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where u — and 1 represent the first and sec- 
ond Cu-O layers, respectively, a is the spin index of 
electrons, c u , r<y = -5= J2k c «,k CT e lk r , TV is the num- 
ber of sites in the lattice and is taken as 400 x 
400 in our calculations below, e k = ti(cosk x + 
cosfc y )/2 + <2 cos k x cos k y + t3(cos2k x + cos2k y )/2 + 
t4(cos 2k x cos fcy + cos k x cos 2k y ) /2 + t^ cos 2k x cos 2k y — /i 
with ti_ 5 =-0.5951, 0.1636, -0.0519, -0.1117, 0.0510 (eV). 
These band parameters are taken from APRES mea- 
surement for Bi2212 [14], and the lattice constant a is 
set as o=l. /i is the chemical potential to be deter- 
mined by doping p. The superconducting order param- 
eter A k = Ao(cosfca; — cosk y )/2. V s is the strength of 
nonmagnetic potential. 

In order to solve the Hamiltonian (1), we first 
take the Bogoliubov transformation: c Ujk | = 

X/i/=0,l( — 1) Cki/V'u.kKi c u,— kj. ~ 2i/=0,l 6knV , t*,kf+l) 

where & = |[1 + (-1)"J£], = and 

_E k = ^/e k + A k , then the Hamiltonian (1) can be 
re-expressed in terms of the fermionic creation and 
annihilation operators tp^ kl/ and VVk^- 
We define two-point Green's functions 

<£ V = -F< T t ^um(t)^ x A0) >, (2) 

where T4>(t) denote the Fourier transform of </>(t) in 
Matsubara frequencies, r is the imaginary time. Fol- 
lowing the standard approach [15], the Green's functions 

Gu.kl" ( iu >n) satisfy 
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where the bare Green's functions G k „ = [iu n — 
(— l)".Ek] _1 - From Eq. (3), we have the Green's func- 
tions with u = 
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which are related to the Green's functions with u = 1. 
Obviously, when the tunneling strength i z = 0, Gq kl/ " 
become the bare Green's functions G kl/ . Therefore, any 
impurity (mid-gap) state in the first Cu-0 layer is pro- 
duced through the inter-layer tunneling. In other words, 
the existence of resonant peak in the surface LDOS pro- 
vides a direct evidence of the tunneling splitting in bilayer 
cuprates. Substituting the expression of G^^ v v (iu> n ) 
into Eq. (3) with u = 1, we have two linearly depen- 
dent equations in terms of G" k ^ v and G" k ^ v , i.e. 
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Solving Eq. (5), we obtain the Green's functions with 
u = 1 
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where D k „ = (a k „a k „+i - & k ,yfe ky+ i) _1 and = 

a ki/+iAp' k ^ - bki/A^^^. We note that a u ', k v and 

/3u',k'i/' in the expression of G™^ v v are undetermined. 
Multiplying Eq. (7) by (-l)"£ki/ and £k,+i, respec- 
tively, and then summating k and v, we also have two 
linearly dependent equations in terms of ov.kv and 
Pu'fc'v' (iw n ), which have solutions 
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FIG. 1: (Color online) The LDOS p(r,u) at the sites (0,0), 
(1,0), and (1,1) with different tunneling strength t z . The 
energy gap Ao = 44 meV and the doping p ~ 15%. 
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FIG. 2: (Color online) The critical value t zc of the tunneling 
strength as functions of the energy gap Ao and the doping p. 



Up to now, we have obtained the Green's functions 
G 



ukl" (i^n)- Obviously, the Green's functions have 

kV 



G^y(iiw„), which 



the symmetries G u ^ (±ioj n ) 
enable us to simplify the calculation of the LDOS below. 

The LDOS at the first Cu-0 layer (i.e. u = 0) can be 



expressed as [16,17] 
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Combining Eqs. (4), (7) and (10), we can calculate the 
curves of the LDOS with various physical parameters at 
different sites and images of the LDOS at different bias 
voltages. In Fig. 1, we present p(r,w) at the sites (0,0), 
(1,0), and (1,1) when the tunneling strength t z = 15 meV 
and 50 meV, respectively. When t z — 15 meV, p(r, uS) 
at these sites are identical below energy gap. When t z = 
50 meV, p(r,u>) has a resonant peak at (0,0) and (1,1) 
with u> = —0.8 meV and at (1,0) with u = 0.8 meV. 
The height of the resonant peaks oscillates decreasingly 
with the distance from the site (0,0) just above the Zn 
impurity. We note that these resonant peaks become 
higher with increasing t z . Undoubtedly, these impurity 
resonances are unique features of the tunneling between 
two Cu-0 layers. 

It is out of question that there exists a critical value 
t zc of the tunneling strength t z , at which a resonant peak 
near the Fermi surface is first produced at (0,0). In Fig. 
2, we show the critical value t zc as functions of energy 
gap Ao and doping p. If p is fixed, t zc decreases with 
increasing Ao. We note that when p ~ 20% t zc does 
not vary with A (> 35 meV). If Ao (<75 meV) is fixed, 
t zc decreases with increasing p. For the other band struc- 
tures, including that with the nearest-neighbor and next- 
nearest-neighbor hopping, t zc has similar relations with 
p and A . 

Now we discuss the STM experiments on a bilayer 
cuprate, e.g. Bi2212. It is known that the surface of 
Bi2212 is a Bi-0 layer and these Bi atoms locate just 
above Cu atoms in the first Cu-0 layer. Therefore, the 
electrical current from the STM tip is blocked by Bi atom 
and is approximately proportional to the summation of 
the LDOS at four neraest neighboring sites of the Cu 
atom beneath the STM tip [12]. The blocking LDOS 
(BLDOS) p S TM(r,w) » T,s=±M r + + p(r + 

(5i y , u>)], where i x and i y are the unit vectors along x and 
y directions, respectively. In Fig. 3, we plot the curves 
of psTM(r,w) at different sites. When t z = 15 meV, 
PSTM(r,w) has no resonant peak near the Fermi surface. 
The curves of psTM(r, w) and p(r, ui) at the same site are 
similar below energy gap [see Fig. 2(a)-(c) and Fig. 3(a)- 
(c)]. When t z — 50 meV, psTM(r, w) has a resonant peak 
at (0,0) and (1,1) with u = 0.8 meV and at (1,0) with 
uj = —0.8 meV. Obviously, the highest resonant peak in 
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FIG. 3: (Color online) The predicted BLDOS psTM(r,cj) 
at the sites (0,0), (1,0), and (1,1) with different tunneling 
strength t z . The energy gap Ao = 44 meV and the doping 
p ~ 15%. 
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FIG. 4: The patterns of pstm(i", uj) near the site just above a 
Zn impurity at different energy uj. The radius of a solid circle 
is proportional to the magnitude of pstm(i",w). The energy 
gap Ao = 44 meV, the doping p ~ 15%, and the tunneling 
strength t z = 50 meV. 



Pstm(i", uj) is at (1,0) while that in p(r, uj) is at (0,0). We 
note that the resonant peaks in pstm{^,uj) and p(r, uj) 
at the same site locate at opposite bias voltages (see Fig. 
l(d)-(f) and Fig. 3(d)-(fj). In Fig. 4, we plot the images 
of pstm(?,uj) near the site (0,0) with t z = 50 meV at 
uj = ±0.8 meV. Obviously, ,ostm(i",w) oscillates with the 
distance from (0,0) and has a maximum at (±1,0) and 
(0, ±1). These patterns also have a four- fold symmetry, 
but are opposite to those induced by a Zn impurity in 
the first Cu-0 layer of Bi2212, where the resonant peak 



is highest at (0,0) [11]. 

It is known that the cuprate superconductors are usu- 
ally inhomogeneous [18]. Different regions have differ- 
ent dopings and energy gaps. For Bi2212, the region 
with larger doping has a smaller gap [19]. From Fig. 
2, we can see that the critical tunneling strength t zc to 
produce a resonant peak through the interlayer tunnel- 
ing depends strongly on doping and energy gap. The 
smaller t zc , the easier the resonant peak to be observed 
by the STM experiments. In order to sec more clearly the 
tunncling-mcdiated resonant peak near the Fermi sur- 
face, we should choose the regions of Bi2212 with larger 
doping. This is consistent with the choice of samples 
in ARPES experiments [9,10]. In other regions with the 
tunneling strength smaller than the critical value, there is 
no resonant peak in the LDOS. From viewpoint of statis- 
tics, the numbers of the Zn impurities in the first and 
second Cu-0 layers in Bi2212 are equal. Therefore, the 
number of the patterns of pstm(i", uj) in Fig. 4 is much 
less than that induced by the Zn impurities in the first 
Cu-0 layer in Bi2212 [20]. 

In conclusion, we have solved analytically the Hamilto- 
nian (1) by the Green's function method and calculated 
numerically the surface LDOS with different physical pa- 
rameters. We have found the critical tunneling strength, 
above which the impurity resonant peak appears near 
the Fermi surface. The LDOS pattern is opposite to that 
produced by the Zn impurities in the first Cu-0 layer of 
Bi2212 at the resonant bias voltage. The sign of the bias 
voltage at the resonant peak alternates with the distance 
from the site just above the Zn impurity in the second 
Cu-0 layer. We hope that the tunneling-mediated impu- 
rity resonant peak and the corresponding LDOS pattern 
could be verified by the future STM experiments. 

The authors would like to thank Professor S. H. Pan 
for useful discussions. This work was supported by the 
Texas Center for Superconductivity at the University of 
Houston and by the Robert A. Welch Foundation under 
the Grant no. E-1411. 



[1] C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

[2] A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 

Phys. 75, 473 (2003). 
[3] S. V. Kivelson et al, Rev. Mod. Phys. 75, 1201 (2003). 
[4] E. Demler, W. Hanke, and S.-C. Zhang, Rev. Mod. Phys. 

76, 909 (2004). 
[5] P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 

78, 17 (2006). 

[6] A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. 

Phys. 78, 373 (2006). 
[7] 0. Fisher et al, Rev. Mod. Phys. 79, 353 (2007). 
[8] O. K. Anderson et al, J. Phys. Chem. Solids 12, 1573 

(1995). 

[9] D. L. Feng et al, Phys. Rev. Lett. 86, 5550 (2001). 



[10] Y.-D. Chuang et ai, Phys. Rev. Lett. 87, 117002 (2001). 
[11] S. H. Pan et ai, Nature 403, 746 (2000). 
[12] J.-X. Zhu, C. S. Ting, and C.-R Hu, Phys. Rev. B 62, 
6027 (2000). 

[13] I. Martin, A. V. Balatsky, and J. Zaanen, Phys. Rev. 

Lett. 88, 097003 (2002). 
[14] M. R. Norman, M. Randeria, H. Ding, and J. C. Cam- 

puzano, Phys. Rev. B 52, 615 (1995). 
[15] G. D. Mahan, Many-Particle Physics (New York: 



Kluwer, 2000). 
[16] Degang Zhang and C. S. Ting, Phys. Rev. B 

100506(R) (2003); Phys. Rev. B 69, 012501 (2004). 
[17] Degang Zhang, New J. Phys. 9, 256 (2007). 
[18] S. H. Pan et ai, Nature 413, 282 (2001). 
[19] J. E. Hoffman et ai, Science 297, 1148 (2002). 
[20] S. H. Pan, private communication. 



